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ABSTRACT 

We study the dynamical evolution of a test particle that orbits a star in the presence of an exterior 
massive planet, considering octupole-order secular interactions. In the standard Kozai mechanism 
(SKM), the planet's orbit is circular, and so the particle conserves vertical angular momentum. As 
a result, the particle's orbit oscillates periodically, exchanging eccentricity for inclination. However, 
when the planet's orbit is eccentric, the particle's vertical angular momentum varies and its Kozai 
oscillations are modulated on longer timescales — we call this the eccentric Kozai mechanism (EKM). 
The EKM can lead to behavior that is dramatically different from the SKM. In particular, the par- 
ticle's orbit can flip from prograde to retrograde and back again, and it can reach arbitrarily high 
eccentricities given enough time. We map out the conditions under which this dramatic behavior (flip- 
ping and extreme eccentricities) occurs, and show that when the planet's eccentricity is sufficiently 
high, it occurs quite generically. For example, when the planet's eccentricity exceeds a few percent 
of the ratio of semi- major axes (outer to inner), around half of randomly oriented test particle orbits 
will flip and reach extreme eccentricities. The SKM has often been invoked for bringing pairs of 
astronomical bodies (star-star, planet-star, compact-object pairs) close together. Including the effect 
of the EKM will enhance the rate at which such matchmaking occurs. 



1. INTRODUCTION 



Naoz et al. (2011a I recently showed that in a system 
with two planets orbiting a star, the inner planet's or- 
bit can flip from prograde to retrograde and back, and 
can also reach extremely high eccentricities. Therefore, 
starting from a system with two prograde planets on dis- 
tant orbits from their star, the inner planet can both 
flip its orbital orientation and reach high enough ec- 
centricities to be tidally captured by the star. That 
might explain the origin of the ~ 25% of hot Jupiters 
whose orbits are retrograde with respect to the spin of 
their star, as in ferred from Rossiter -McLaughlin mea- 
surements (e.g., |Triaud et al.||2010j). However, while 

is in 



the behavior displayed in |INaoz et al. 



triguing, it has not yet been understood. Furthermore, 



(2011a 



Naoz et al. (2011a) choose initial eccentricities and in- 
clinations that are fairly extreme for planetary systems 
(e ~ 0.6, i ~ 70°), and hence it is not yet clear how 
generic their results are. The goal of this paper is to 
understand and map out the mechanism in a simplified 
system where the inner planet is treated as a massless 
test particle. 

There has been much work in the litera t ure o n the 
orbital dynamics of similar systems. |Kozai| ( 1962 ) con- 
sidered the evolution of an asteroid perturbed by a cir- 
cular Jupiter. He focused on secular interactions, mean- 
ing that interactions are averaged over the orbital phases 
of the asteroid and Jupiter. Kozai found a remarkable 
result: if the asteroid's orbit is sufficiently inclined (be- 
tween 39° and 141°), then it cannot remain on a circular 
orbit. Instead, its eccentricity and inclination oscillate 
periodically. Furthermore, if the asteroid is highly in- 
clined (~ 90°), then its eccentricity will grow over the 
course of a single Kozai oscillation to ~ 1. 

The Kozai mechanism has been applied to a variety of 
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astronomical systems. For example, in triple star sys- 
tems, if the outer star is inclined it can force Kozai oscil- 
lations in the inner binary, increasing the inner binary's 
eccentricity until ti dal dissipation circularizes it into a 
very tight orbit (e.g . , Egglcton & Kiseleva-Eggleton 2001 



Fabrycky fc Tremaine|2U07| . Similarly, it a Jupiter-mass 



planet forms around a star with an inclined compan- 
ion star, that companion can force Kozai oscillations in 
the planet u ntil tides damp its orbit, forming it into a 
hot J upiter ( |Wu fc Murray|)2003j |Fabrycky fc Tremainej 



2007 1 . Analogous scenarios have been pro posed for merg- 



ing supermassive black holes (Bla es et al. 2002 1 and 
stellar-mass compa ct objects (Thompson 2010). See also 



Naoz et al. (2011b I for a review of other applications of 
the Kozai mechanism. 

In the case considered by Kozai, the asteroid's ver- 
tical angular momentum (=const x\/l — e 2 cos i) is con- 
served because Jupiter's orbit is circular. Since this is 
a two degree of freedom system (eccentricity and incli- 
nation) with two conserved quantities (vertical angular 
momentum and secular energy), it is integrable. But 
when Jupiter's orbit is eccentric, the dynamics can no 
longer be solved analytically. As Kozai noted, "[wjithout 
the aid of a high-speed computer, it is rather difficult to 
estimate the effects of Jupiter's eccentricity." 

The effect of an eccentric perturber on Kozai oscilla- 
tions has been con sidered in a number of example cases. 



Harrington ( 1969 1 integrates the secular octupole equa- 



tions tor tour triple star systems, and shows that in one 
case there is a new resonance; in a second case there 
is chaos; and in the other two cas es no new interesting 
effects arise. Holman et al. (1997) perform a direct nu- 
merical integration tor a planet forced by a star, but with 
their parameters (e ~ 0.008 in our notation, see below), 
the perturber's finite eccentricity only leads to a narrow 
zone o f chaos around the Kozai separatrix. |Ford et al.| 
( 2000 1 integrate the secular octupole equations for a sys- 
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tern of two eccentric planets, where the outer planet has 
eccentricity = 0.9, and find that the inner planet's Kozai 
oscillations are modul ated on long t imesc ales, leading to 
extreme eccentricities 
examples for t riples of 



pool 



examine some 



Blaes et al. 

supermassive black holes. [Naoz| 
et al. ( 2011a|b ) examine some examples for two planets. 



Yet a systematic exploration of the effect of an eccentric 
Kozai perturber is lacking. 

In a number of studies the effect of an eccentric per- 
turber on Kozai oscillations has be en neglected, even in 
situations where it should not be (Naoz et al. 2011b). 
Part of the reason is that treatments of Kozai oscilla- 
tions often expand the secular Hamiltonian to leading 
order in the ratio of semimajor axes, i.e. to quadrupole 
order. At that order, the exterior body's argument of pe- 
riapse does not appear, which implies that the particle's 
vertical angular momentum is conserved^] To see the ef- 
fect of the eccentricity, one must work to higher order 
in the ratio of semimajor axes, i.e. at least to octupole 
order. 

In this paper, we extend Kozai's work to the case of 
an eccentric planet and map out the resulting behavior, 
which we call the eccentric Kozai mechanism (EKM). 
To do so, we use the secular octupole Hamiltonian of 
|Ford et ah] (|200 0"|), which was ap parently first der ived 
by [Harrington ( 1968]) ( as cited by Harrin gton 19691 see 

Me 



and |Krymolowski &: IVlazeh 1999 for 



also |Marchal||199 

derivations) . 

After the work described here was completed, we learnt 
of a similar paper being prepared by Katz, Dong, and 
Malhotra. 

This paper is organized as follows. In Section [2] we 
present the equations of motion, relegating their deriva- 
tion to the Appendix. In Section [3j we review the SKM. 
In Section |4j the heart of this paper, we map out the 
EKM. We summarize in Section O 



2. EQUATIONS OF MOTION 

We solve for the orbit of a massless test particle in the 
presence of an exterior massive planet, including only 
secular interactions expanded to octupole order. The 
planet is on a fixed eccentric orbit, and the particle's 
orbit is specified by four variables, 



{e,w,9 = cosi,£!} 



(1) 



which are its eccentricity, argument of periapse, inclina- 
tion (or its cosine), and longitude of ascending node rel- 



ative to the planet's periapse (e.g., [Murray fc Dermott 
20001. In the Appendix, we use the secular octupole 



Hamiltonian that has been published in the literature to 
derive the particle's equations of motion. As we show in 
the Appendix, although that published Hamiltonian has 
had its nodes eliminated, one can still use it to derive the 
full test particle equations of motion — even the equation 
that requires the nodes. 

We summarize the equations here. Defining the parti- 
cle's (scaled) total angular momentum and vertical an- 

2 The fact that the exterior body's argument of periapse docs 
not appear to quadrupole order has been called " a happy coin- 



cidence" because it makes t he system integrable jLidov fc Ziglin] 
|1976| |Laskar fc Boue]|2010} . However, it is perhaps more ot an 
unhappy coincidence m view of the fact that it has misled some 
researchers into neglecting the role of the planet's eccentricity. 



gular momentum a^] 



J= y/T- 
J z = 9^/l^ 



,2 



(2) 
(3) 

the equations of motion may be expressed as partial 
derivatives of an energy function F(e, uj, 9, O) via 
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where t is proportional to time (Eq. | A8| ) . These are 
Hamilton's equations for the two pairs ol canonically con- 
jugate variables {J, uj} and {J Z ,Q}, except that we ex- 
press F as a function of the non-canonical variables e and 
9 — for that reason, we call F the energy function rather 
than the Hamiltonian. Of course, F is a constant of the 
motion. It is given by a quadrupole and an octupole term 

F = F qu + eF oc , (8) 



where the constant 



(a/opi)e p i 



1 



(9) 



-pi 



here a/a p i is the ratio of semimajor axes (inner to outer) 
and e p i is the planet's eccentricity. The quadrupole piece 
is 

^ qu ee + 9 2 + |e 2 2 + \e\\ 9 2 ) cos(2 W ) , (10) 

after dropping an irrelevant constant, and the octupole 
term is 



F n , 
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)9' 2 - 159 3 ) cos(cj + ft) 
) 2 + 159 3 ) cos(w - fi) 
2 + 9 3 ) cos(3w - 0) 
- 9 3 ) cos(3w + fi)l (11) 



The only adjustable parameter in the equations of mo- 
tion other than the initial conditions is the constant e. 
That constant encodes the properties of the planet]^] In 
the standard Kozai mechanism (SKM), the planet's or- 
bit is circular (e = 0). Hence F is independent of f2, and 
thus J z is a constant of the motion. But if the planet's 
eccentricity is not zero (e > 0), then the octupole term 
allows J z to change. If the planet is either nearly circular 
(e p i <C 1) or distant (a/a p \ <C 1), then e is very small, 

3 Our variables are related to the Delaunay variables {G, g, H, h} 
that ar e often u sed in treatments of K ozai oscil lations via J = 
G/( m^/G NM*a), u = g, J z = H/Oms/UffM^a), U = h (Eqs. 
|A3| -[A4|). 

The planet's properties also ente r in a trivial way through the 
scaling between t and time (Eq. IA81). 
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Fig. 1. — SKM (e = 0): The three left panels show trajectories 
with J z = 0.2 and various values of F. Trajectories with the same 
Jz overlap in the e-6 plane (bottom panel); we have slightly offset 
these for clarity. The right panels show the same, but with J z = 
0.6. For circulating trajectories, the minimum e and 9 occur at u = 
0. The dashed hor izontal lines show the cosine of the critical Kozai 
inclination, y3/5. The separatrix — which separates circulating 

from librating trajectories — always has eo = and 1 9 w /2 1 = \/3/5 



and the evolution is typically similar to the SKM (except 
at extremely high inclinations). But as we shall show, if 
e > 0.01, then the octupole term can lead to qualitatively 
new behavior that is not found in the e = limit. 

3. THE STANDARD KOZAI MECHANISM (e = 0) 

We review the SKM to set the stage for the EKM. 
For the SKM, the planet's orbit is circular (e = 0, i.e. 
F = F qxl ). Hence there are two constants of the motion, 
F and J z , and the motion is regular. Each trajectory may 
be labelled by the values of F and J z . Inserting Equation 
(13| into Equation ( |10| immediately determines e as a 
function of a;, as well as 8 as a function of u>. The three 
left panels of Figure [T] show sample trajectories, with 
fixed J z = 0.2 and various values of F; the right panels 
show the same, but with J z — 0.4. For our purposes, the 
most important properties of the SKM are as follows: 



• Because J z =constant, each trajectory traces out 
a curve in the e-6 plane, which we call a "Kozai 
curve." 

• There are two classes of trajectories, librating and 
circulating. On circulating trajectories, e and \6\ 
are smallest at Li = 0, and they are largest at u> = 
±7r/2. The separatrix has eo = el^o = fj^j and 

\6±*/2\ = 

Throught this paper, we denote values at u = with subscript 



-i — I 1 1 — r- 




r 



Fig. 2. — SKM (e = 0): The left panel shows trajectories with 
fixed F and various values of J z ("Kozai curves"). The left-most 
tips of the Kozai curves (with co-ordinates {eo,6*o}) trace out 
another curve, which we call the "energy curve." That energy 
curve is shown in green in the right panel; it is the level curve of 
^qu.o = ^qu|u=o = 0.16. The right panel also shows some other 
energy curves. These intersect the 6>o axis at -F qu ,o and the eo 
axis at ^/F , qu ,o/2 

• On a trajectory that has eo <C 1, the largest e is 

5 



S ±7T/2 



i - -el 



1 - i;J'i , 
3 z 



to leading order in e\ , when | 6q \ 



(12) 

(The cor- 
responding inclinations are the critical Kozai angles 
cos' 1 ±^3/5 = 39° and 141°.) Therefore when 
1 6*o I "C 1 (i.e. inclination close to 90°), the largest 
eccentricity is nearly unity. 

• Given F and J z , the minimal e on a circulating 
trajectory satisfies 



4 = \(F- 



(13) 



Figure [2j left panel, shows a sequence of Kozai curves 
with the same F and differing J z . The left boundaries 
of the sequence (i.e., the values {eo,#o} for each Kozai 
curve) trace out a curve in the e-6 plane, which we call 
an "energy curve." That curve is plotted in green in the 
right panel of Figure [2] An energy curve is a curv e of 
constant F qu | u= o. Its form is given by Equation_ ( 13 1 
with Jl = e\(\ — 6*q). The right panel of Figure 
shows other energy curves with different energies 



also 



qu- 



4. THE ECCENTRIC KOZAI MECHANISM (e > 0) 

When the planet is eccentric (e > 0) there is only a sin- 
gle conserved quantity, the secular energy F . Therefore 
the particle's trajectories are more complicated, and can 
even be chaotic. Figure [3] shows two sample trajectories 
in the e-6 plane for e = 0.01. Both trajectories have the 
same energy F = 0.16, but different initial conditions 
consistent with that energy. Rather than being confined 
to a single Kozai curve, the particle evolves from one 
Kozai curve to another. As long as e < 1, all of these 
Kozai curves have nearly the same F qu , and therefore the 
trajectories in Figure [3] follow along the tracks displayed 
in the left panel of Figure |2l 

Trajectory 'a' in Figure 13] evolves through 6 = 0. Its 
orbit flips from prograde (# > 0) to retrograde (6 < 0) 
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Fig. 3. — EKM: Two sample trajectories ('a' and 'b') with F = 
0.16, e = 0.01, and different initial conditions. When the planet 
is eccentric, the test particle no longer traces out a single Kozai 
curve, but evolves from one Kozai curve to another, approximately 
tracing out a sequence of Kozai curves that have a fixed quadrupole 
energy (Fig. [2] left panel). Orbit 'a' flips repeatedly (i.e. crosses 
6 = 0) and reaches extreme eccentricities, and orbit 'b' does not. 
Although only points up to time t = 500 are plotted, longer runs 
yield identical plots, albeit with more densely packed points. In 
particular, orbit 'b' never flips. 

and back again. Furthermore, its eccentricity approaches 
unity. In fact, the flipping of an orbit is closely tied to 
its eccent ricit y reaching unity. This can be seen from 
Equation ( 12 ), which implies that a Kozai curve that has 
9q ~ reaches e ~ 1. (It never reaches exactly e = 1; 
see below). Trajectory 'b' never flips, and its eccentricity 
does not approach unity. 

The left panels of Figurej4]show the temporal evolution 
of 9 and e for trajectory 'a .The bottom left panel shows 
two angles that appear as arguments of cosine terms in 
F oc . When the particle's orbit is prograde, the angle ui + 
fi librates and uj — fl circulates; and when it is retrograde, 
those two angles switch roles. The right panels of Figure 
[4] show the same quantities for trajectory 'b'. Whereas 
trajectory 'a' is regular, 'b' is chaotic. We demonstrate 
that more explicitly below. 

The left panels of Figure [5] are a zoom-in of the 'a' tra- 
jectory at early times. The fast oscillations in 9 and e 
are primarily governed by the quadrupole evolution, as e 
and 9 trace out Kozai curves with F — 0.16 and various 
values of J z . Over the course of a single oscillation, w 
increases from to n or from — ir to . Hence the or- 
bit is always circulating (see Fig. 111. The top left panel 
also shows J z . Whereas in the SKM J z =const., here J z 
changes in a nearly step-wise fashion. There are sharp 
jumps in J z whenever e and 9 change rapidly; these are 
forced by the octupolar contribution F oc . The long-term 
evolution of J z controls the envelopes of both 9 (Fig. [5j 
top-left panel) and e (Fig. [51 middle-left panel; see eq. 
|12| ). Successive maxima of e occur in discrete steps. 
Therefore even when J z crosses through zero, the max- 
imum e is never precisely unity. Nonetheless, as time 
evolves, the maximum e reached approaches closer and 
closer to unity. 

The right panels of Figure [B] show the corresponding 
zoom-in for trajectory 'b'. Thebottom right panel shows 
that w switches from circulation to libration and back 
again. That explains why t he 'b' trajectory is c haotic. 



Similar behavior is shown in Holman et al. ( 1997|) 
A global view of the dynamics is provided by sur 
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Fig. 4. — EKM: Temporal evolution of the two trajectories de- 
picted in Figure|3] In the bottom panels, the angles are cumulative. 
Trajectory 'a' reaches extreme eccentricities when it flips. Its evo- 
lution is regular, with uj + Q librating when it is prograde, and 
u) — f2 librating when it is retrograde. Trajectory 'b' is chaotic, and 
never flips. 
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Fig. 5. — EKM: Zoomed in evolution of trajectories 'a' and 'b' at 
early times (Fig. [3}. In the top panels the re d lin es show that J z 
controls the envelope of both 8 and e (via eq. |12| ). In the bottom 
panels, the angles are not cumulative. 
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Fig. 6. — EKM, e = 0.01. Surfaces of section with various values 
of F; points are plotted whenever uj = 2nn for integer n. For each 
F, if the values of 8q were plotted a gai nst the value of eo, they 
would lie along an energy curve (Fig. [2] right panel). The colors 
in that figure correspond to the colors in this one. 



section. Figure [6] maps out the behavior when e = 0.01. 
We plot a point whenever u = 2nir for integer n This 
may be interpreted as follows. Consider trajectory 'a' of 
Figure [3] When to = 2mr that trajectory hits the green 
energy curve depicted in the right panel of Figure [2j At 
those times we plot in Figure [6] the value of e = eo versus 
Q, and 9 = 8q versus f2 (curves labelled 'a', middle pan- 
els). Also shown in those green middle panels are other 
trajectories with the same energy F — 0.16, including 
trajectory 'b'. If all of the green points were plotted 
against each other in the eo-#o plane, they would trace 
out the energy curve labelled 0.16 in Figure [2] Equiva- 
lently, one can imagine that that energy curve is extended 
out of the plane of the paper into a half-cyliner, with the 
third dimension being the value of ft. The F = 0.16 
surfaces of section of Figure [6] cover the surface of that 
half-cylinder. Of course, for a given energy the maxi- 
mum | 6q I is yF, and the maximum eo is \J F/2. Beyond 
those values, the energy curve does not exist (Fig. [2] 
caption). 

The other panels in Figure 6] may be interpreted sim- 
ilarly, with each pair of panels at fixed F lying along a 
single energy curve of Figure 2] (with corresponding col- 
ors). We therefore now have a virtually complete view of 
the e = 0.01 dynamics. 

Perhaps the most dramatic new behavior caused by 
the octupole term is that orbits can flip orientation (i.e. 
cross (9 = 0), and, as a consequence, reach arbitrarily high 
values of eccentricity. From Figure |6j it is immediately 
clear which orbits exhibit this behavior, and which do 
not, at e = 0.01. For example, orbits with eo = will 

6 Trajectories only appear on our sections while u> is circulating, 
because ui is never equal to 2nn during librations (Fig. 111. 
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Fig. 7. — First Flipping Orbit. Each curve marks the first flipping 
orbit that can occur for a given value of e. Above each curve (i.e., 
at low inclination) none of the orbits flip, and below it some do. 

always flip provided that \8q\ < 0.2 (i.e., io > cos -1 0.2 ~ 
80° for prograde orbits). For eo ^ 0, we infer from the 
blue (F — 0.25) panels that even orbits with 9o as large 
as 0.4 (io as small as 66°) can flip. 

A second new behavior is the appearance of chaos. 
Regular orbits appear as curves in the surfaces of section, 
while chaotic orbits appear as a smattering of points. For 
example, it is apparent from Figure [6] that trajectory 'a' 
is regular and 'b' is chaotic. Chaos always occurs near 
eo = 0. That is because the Kozai separatrix always has 
eo = (Section [3]), and chaos is caused the crossing of 
the Kozai separatrix, when u transitions from librating 
to circulating and back to librating (e.g., Fig.j5l). 

Figure [7] is the main result of this paper. It summa- 
rizes where flipping orbits occur for various values of e. 
Towards the top of the plot (large 9q and small io), there 
are no flipping orbits. But as 9q decreases, an increasing 
number of orbits can flip. Each curve marks, for a given 
e, the location at which the first flipping orbit occurs, if fi 
is judiciously chosen. Above the curve there are no flip- 
ping orbits, and below there are some — typically around 
half of the orbits below the curve can flip. We made 
this plot by running a suite of numerical simulations ini- 
tialized with 9 = and uj = 0, and scanning through 
values of initial e and f2 to determine which flipping or- 
bit reached the largest values of 9g for a given initial e, 
and then plotting the extreme values of eo and 9q for that 
orbit. (Note that all such trajectories initialized with the 
same initial e — eo share the same energy curve.) 

A striking feature of Figure [7] is the fairly sudden tran- 
sition when e exceeds a few percent. For example, focus- 
ing on orbits with eo = 0, when e = 0.01 flipping orbits 
can occur only if \9 Q \ < 0.2 i.e., io > 80° (for prograde), 
as stated above. But for e = 0.03 flips can occur for 
1 6*o I ~ 0-5, i.e. io ~ 60° (for prograde). Note that if test 
particle orbits are randomly oriented, they will be uni- 
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formly distributed in and therefore half of the orbits 
will have 9 < 0.5. Of course, this is only a rough esti- 
mate for the fraction of flipping orbits, since it depends 
also on the values of the other variables (e, ui, f2). 

The curves in Figure [7] rise as eo increases from 0, and 
then suddenly turn over. The reason for the sudden turn 
over is apparent from a glance at the surfaces of sec- 
tion (Fig. 6|. In the blue (F — 0.25) section, there are 
flipping orbits centered both on fl = and on ft = tt, 
with the latter ones reaching to a larger \6q\. But in the 
magenta (F = 0.36) section, there are no more flipping 
orbits centered at Cl = tt. It is that discontinuity that 
gives rise to the sharp break in Figure [7] 

5. SUMMARY 

In this paper, we map out the behavior of a test parti- 
cle that is forced by a massive eccentric planet, working 
under the secular octupole approximation. 

In Section[2j we give the evolutionary equations. These 
are derived in the Appendix, where we show how one 
may use the published Hamiltonian — which has had its 
nodes eliminated — to derive the full equations of motion. 
There is a single parameter ,e, that characterizes devi- 
ations from the SKM (Eq. [9]). The SKM is recovered 
when e -C 1, i.e. when the planet is either nearly circular 
or very distant. 

In Section [3j we review the SKM. In the e-9 plane 
(where 9 = cosz), the particle traces out a Kozai curve 
because its J z is fixed. A sequence of Kozai curves with 
the same energy F outline an energy curve, which is a 
curve of constant secular energy, evaluated at u = 0. 

In Section [4j we map out what happens when the 
planet's orbit is eccentric. Our principal results are: 

• In the e-9 plane, a single trajectory evolves from 
one Kozai curve to another, all of which abut the 
same energy curve (Fig. [3]). 

• For trajectories where the orbit flips, the ec- 
centricity becomes extremely high, and as time 
evolves, the maximum eccentricity reached ap- 
proaches closer and closer to unity. That is because 
every time 9 crosses through zero, |0o| becomes 
very small. And given enough zero-crossings, the 
smallest \9 \ reached becomes arbitrarily small, and 
hence the maximum eccentricity becomes arbitrar- 
ily close to unity (as determined by Eq. 12 ). 



to circ ulating, and back. See also |Holman et al. 
(p97l). 



• The vertical angular momentum J z varies gradu- 
ally in the EKM, unlike in the SKM where it is 
constant. Its evolution controls that of both 9 and 
e (Fig. [4]). Even when the value of e is small, 
and hence the temporal evolution of J z is slow, if 
J z crosses through zero, the dramatic flipping and 
extreme eccentricity behavior occurs (Fig. [31 left 
panel) . 

• Chaotic behavior occurs when the Kozai separatrix 
is crossed. In that case, to transitions from librating 



• Surfaces of section, taken whenever lu = 2mr, pro- 
vide an graphical global map of the dynamics. 
Each energy curve gives rise to two panels in Fig- 
ure [6] (top and bottom). By stepping through the 
energy curves, one acquires a complete map of the 
dynamics for a given value of e. From the sur- 
faces of section, one can immediately see which or- 
bits exhibit flipping and extreme eccentricity, and 
which do not. One can also see which orbits are 
chaotic and which are not. All four combinations 
occur: flipping and non-chaotic ('a'), flipping and 
chaotic, non-flipping and non-chaotic, non-flipping 
and chaotic ('b'). 

• Figure [7] is the main result of this paper. For each 
e, it gives the curve in the eo-#o plane at which the 
first flipping orbit occurs (along with the accompa- 
nying extreme-eccentricity behavior). Above the 
curve no flipping orbits occur, and below it an in- 
creasing number occur. There is a fairly sharp tran- 
sition in the smallest inclination at which orbits 
can flip. For e ~ 0.01, the inclination must exceed 
~ 80° to flip (for prograde orbits with eo = 0); but 
when e ~ 0.03, it need only exceed ~ 60°, implying 
that around half of randomly oriented test particle 
orbis will flip in that case. 

We have made a number of approximations in order 
to simplify the analysis. First, we work under the sec- 
ular octupole approximation, which is poor when e is 
too large because in that case the particle comes close 
to the planet. We have made a few comparisons with 
direct integra tions, some of which are presented in |Naoz| 
et al. (2011b). Thus far the agreement is quite good, but 
more should be done to determine the domain of valid- 
ity of the secular octupole approximation. Second, we 
take the inner orbiting body to be massless. While that 
is appropriate for a planet being forced by a companion 
star, it is not for t wo planets, which i s the main motiva- 
tion for this paper ( Naoz et al.|2011a ). We are currently 
extending the work presented here to two massive plan- 
ets; in that case, there is a gen eralized e that includes 
the masses (Naoz et al. 2011b). And third, we ignore 
all physical effects beyond Newtonian gravity, including 
tides and general relativistic precession. These will play 
a role when the inner planet's eccentricity approaches 
unity. 
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APPENDIX 

OCTUPOLE SECULAR EQUATIONS OF MOTION FOR A TEST PARTICLE 

The secular interaction energy between an interior and exterior planet (unprimed and primed, respectively) can be 
written as 

TTWd! 

E eon = -G N — — a 2 f(a,e',e,g',g,9) , (Al) 
a 

where 9 = cos / and the other quantities are the mutual inclination /, argument of periapse g, eccentricity e, semimajor 
axis a, ratio of semimajor axes (inner to outer) a, mass m, and Newton's constant G jv- (The notation here differs 
slightly from the body of the paper in order to connect to previous treatments of the octupole. See footnote [3] for the 
relation between the two notations.) The subscript "eon" refers to the elimination of the nodes, i.e. the longitudes of 
ascending nodes (h, h!) have been chosen to satisfy h — h' = tt which is why they do not appear in the energy. 

We seek the equations of motion for the inner planet in the limit that m — > 0. The equations of motion can be 
expressed as Hamilton's equations for the canonical Delaunay variables {g, G; h,H} where 

G = m^/G N M*ay/l - e 2 , H = Gcosi, (A2) 

with M* the mass of the central star. Note that in secular theory, a is constant; note also that in the test particle 
limit, i = I. 

Setting the Hamiltonian H = _ E eon will give the co rrect equations of motion of g, G, and h, but not for H because 



the nodes have been eliminated (Naoz et al. 2011b). We shall prove that the correct procedure for obtaining the 
Hamiltonian is simply to replace g' — > tt — h in the above expression for E con , as well as e' =const., e 2 — > 1 — 
G 2 / (m 2 G nMqo), and 8 — > H/G. The relation g' — > tt — h might appear surprising because one might have expected 
that g' =const when the inner particle is massless. But in truth h' is undefined because the reference plane is aligned 
with the outer orbit. Therefore, the outer planet must only have g' + /i'=const, and we may choose without loss of 
generality the constant to equal zero. Hence, elimination of the nodes implies g' — n — h 

Proof 

Our proof is based on the triangle equality G 2 + G' 2 + 2GG'9 = (H + H') 2 , which expresses the square of the vertical 
angular momentum in two equivalent ways. Taking the time derivative of the triangle equality, and noting that both 
G and G' satisfy Hamilton's equations, i.e. dG/dt — —dE con /dg and dG'/dt = — dE con /dg' , we find that to leading 
order in m, dG'/dt — —^G9, i.e. dH/dt = dE con /dg' which proves that if we replace g' — > n — h in E con , we will end 
up with the correct Hamilton equation for H. 

Equations of Motion 

We may re-scale all momenta by an arbitrary constant without changing the equations of motion as long as we 
rescale the Hamiltonian by the same constant. Therefore defining 

G 



J= = y 7 ! - e 2 (A3) 

my/GNM^a 

J z = JcosI = 9\/l - e 2 , (A4) 

the rescaled Hamiltonian is 

H(J,g;J z ,h)= (A5) 

= -^-n,a 3 f(a,e\e,g\g,9) , (A6) 
where ft* is the orbital angular speed of the test particle around the star. 



The expression for / (defined via eq. IAl| ) is derived, e.g., in Ford et al. (2000) to octupole order. We re-express 
their Equation (22) in the limit m, m' <C 



f ~ 8(l-e' 2 ) 3 / 2(Fqu + l-e'2 Foc) ' 



(A7) 
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and [TT1 with g = uj and h = f2) after correcting the sign of the 
7T — h 



where F qu and F oc are displaye d above (Eqs 
octupole term ( Ford et al.||2004[ ) and replacing g 

The equations of motion are Hamilton's equations for Hamiltonian ( A6). These are displayed explicitly in the body 
of the paper (Eqs. [4]- [7]) after defining the rescaled time 



Tn' 1 3 1 



(A8) 



where T is the true time. 



